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ABSTRACT 

We present a survey for 2 ~ 3 Lyman break galaxies (LBGs) associated with 
damped Lyman a systems (DLAs) with the primary purpose of determining the 
DLA-LBG cross-correlation. This paper describes the acquisition and analysis of 
imaging and spectroscopic data of 9 quasar fields having 11 known 2; ~ 3 DLAs 
covering an area of 465 arcmin^. Using deep m'BVRI images, 796 LEG candidates 
to an apparent R^s magnitude of 25.5 were photometrically selected from 17,343 
sources detected in the field. Spectroscopic observations of 529 LEG candidates 
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using Keck LRIS yielded 339 redshifts. We have conservatively identified 211 
z > 2 objects with {z) = 3.02 ± 0.32. We discuss our method of z ~ 3 LBG 
identification and present a model of the m'BVRI photometric selection function. 
We use the 339 spectra to evaluate our w'BVRI ^ ~ 3 Lyman break photometric 
selection technique. 

Subject headings: galaxies: high redshift — quasars: absorption lines — galaxies: 
photometry — galaxies: formation 

1. INTRODUCTION 

Over the last decade, high redshift surveys to detect Lyman break galaxies (LBGs) 
by photometric pre-selection, (Steidel ct al. 1996a,b; Lowenthal et al. 1997), have been 
successful in compiling catalogs of ~ 1000 galaxy spectra at 2; ~ 3. This has become 
an efficient technique to obtain statistical information of high redshift galaxies via their 
emission properties. Analysis of the LBG spatial distribution shows significant clustering 
which supports the belief that they are biased tracers of mass and form in overdense regions. 
The LBG auto-correlation function (Steidel et al. 1998; Adelberger et al. 1998, 2003) and 
best- fit star formation models (Shapley et al. 2001) suggest LBGs are massive objects that 
most likely evolve into early-type galaxies by the present day. LBG surveys, such as the 
spectroscopic sample, are magnitude limited and only probe the higher end of the proto- 
galactic mass distribution. 

By contrast, quasar (QSO) absorption line systems represent a population of high red- 
shift protogalaxies whose detections are not biased by their intrinsic magnitude. Their 
detection is dependent only on the strengths of their absorption line features and the bright- 
ness of the background QSO. The Damped Lyman a systems (DLAs) are the highest column 
density absorbers, defined to have N(H1) > 2 x lO^'' atoms cm~^ as established by Wolfe 
et al. (1986). DLAs are drawn from the bulk of the proto-galactic mass distribution and 
harbor the majority of the neutral hydrogen present in the universe. The agreement between 
the comoving mass density of neutral gas at ^ > 2 and the mass density of visible stars in 
local galactic disks (Wolfe et al. 1995) lends credibihty to the idea that DLAs evolve to form 
present-day galaxies such as the Milky Way (Kauffmann 1996). However, this agreement in 
baryon mass cannot differentiate between passive models in which DLAs are protogalaxies 
having properties similar to those of local galaxies with constant comoving densities from 
CDM models in which DLAs are numerous, low-mass objects with comoving densities that 
decrease with time. 
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These two interpretations have been substantiated largely from detailed analyses of 
DLA gas kinematics. Prochaska & Wolfe (1997, 1998) and Wolfe & Prochaska (2000a,b) 
best described the asymmetries in the absorption line profiles using large thick centrifugally- 
supported disks drawn from a Schecter luminosity/mass function that obey the Tully-Fisher 
relation. Conversely, hierarchical CDM scenarios confine DLAs to compact low-mass sub-L* 
systems. Models of irregular merging clumps are shown to reproduce the observed absorption 
line kinematics as well (i.e. Haehnelt et al. 1998; Haehnelt et al. 2000). A rehable determina- 
tion of the mass of typical individual DLAs is needed to break the degeneracy between these 
interpretations and provide a vital parameter for numerical models of galaxy formation. 

DLAs offer a wealth of information from their absorption line characteristics but their 
baryonic mass cannot be determined without knowing their typical cross-sectional area. 
The incidence of DLAs along quasar lines-of-sight only determines the product of their 
comoving densities and cross-sections. It is possible, however, to infer the dark-matter 
mass of typical DLAs in the context of CDM cosmology. Hierarchical structure formation 
models predict that the most massive galaxies will preferentially form clustered together 
in regions of mass overdensity. Under the assumption of a Gaussian distribution of density 
fluctuations in the density field, nearly all of the density fiuctuations that give rise to massive 
galaxies (with the exception of those on the high-end tail of the distribution) have initial 
density contrasts below that necessary for collapse against the Hubble flow. However, these 
density fluctuations receive a boost in density when superimposed onto a larger wavelength 
fluctuation from an over-dense region. In this situation, the number of the fluctuations 
exceeding the threshold of collapse that correspond to massive galaxies is increased with 
the tendency of these fluctuations to be clustered near the peak of the underlying overdense 
region. Conversely, density fiuctuations that give rise to low-mass galaxies have a relatively 
high initial density contrast and are predicted to form more uniformly throughout space. 

An enhancement in the clustering of galaxies is referred to as bias (Kaiser 1984). There- 
fore, the measurement of the bias of a population of galaxies provides a means to infer their 
mass distribution. The bias of a particular sample of galaxies can be derived from their 
spatial auto-correlation or by their cross-correlation with another population of galaxies 
with known bias. DLAs are detected in absorption in quasar lines-of-sight and their sparse 
sampling at 2; ~ 3 makes the measurement of the DLA auto-correlation function difficult. 
Since spectral analysis from previous surveys for LBGs in emission has yielded the LEG 
auto-correlation function and estabhshed the LEG galaxy bias at 2; ~ 3 (Adelberger et al. 
2003), an effective way, and perhaps the only way, to determine the DLA bias is through the 
DLA-LEG cross-correlation. 

In general, the Poisson probability that two volume elements, dVi and dV2, in the held 
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separated by a distance r are occupied by a DLA and an LBG is: 



dP{r) = UDLAriLBG dVoLAdVLBG, 



(1) 



where udla and ulbg are the mean number densities of each population. The presence 
of a DLA in one volume element enhances the probability of the presence of an LBG in 
another volume element by an amount described by the DLA-LBG cross-correlation function 
iDLA-LBcir). Therefore, the probability of finding an LBG in a volume element at a distance 
r from a known DLA is expressed as: 



At high redshift, where a linear bias to the underlying dark matter correlation function ^dm 
(i.e. Peacock 1997) is a reasonable model, the DLA-LBG cross-correlation is: 



where b^iA and 6lbg are the bias parameters for the two populations. Since the LBG bias 
at z ~ 3 is known, and once the DLA-LBG cross-correlation is measured, the dark matter 
mass of typical individual DLAs can be inferred directly from the resulting value of the DLA 
bias using the approach of Mo et al. (1998). 

To measure the DLA-LBG cross-correlation function, we use a photometric selection 
technique similar to previous successful surveys, Steidel et al. (1998), to efficiently detect 
LBGs in QSO fields with known DLAs. We are able to construct the DLA-LBG spatial cross- 
correlation function from both angular and spectroscopic measurements. Recently, Bouche & 
Lowenthal (2003) estimated the two-dimensional DLA-LBG angular cross-correlation func- 
tion in three wide fields using photometric redshifts. They found a positive DLA-LBG 
clustering amplitude to a > 95% confidence level based on Monte Carlo simulations. Their 
photometric redshift measurements constrained their results to spatial clustering within red- 
shift slices on the order of their photometric redshift uncertainties (A^ ~ 0.2). Conversely, 
Adelberger et al. (2003) measured the mean overdensity of LBGs within cylindrical cells with 
a height Az = 0.025 and radius of up to = 265" centered on the four DLAs in their spec- 
troscopic sample at 2; ~ 3. They compared this mean overdensity with the mean overdensity 
of LBGs within similar cylindrical cells centered on LBGs drawn from their complete sample 
that had similar redshifts as two of their DLAs. Prom their analysis, they argued that the 
three-dimensional DLA-LBG cross-correlation is weak and concluded that DLAs and LBGs 
do not have similar spatial distributions. Since the statistical accuracy of both investigations 
was limited, neither confirms, nor rules out a significant value of the DLA bias. 

Because DLAs may be a population of protogalaxies having a range of masses, a statis- 
tical determination of the DLA-LBG cross-correlation function from multiple DLA fields is 



dP{r) — nDLAriLBG [1 + CDLA-LBGit^)] dVoLAdVLBG- 



(2) 



ioLA-LBG — boLA blBG ^DM 



(3) 
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required. This paper presents the results of u'BVRl photometric selection of LBG candidates 
in a survey of 9 QSO fields containing 11 known z ~ 3 DLAs. We address the efficiency 
of using m'BVRI filters to detect LBGs and apply 339 spectra to calibrate the m'BVRI pho- 
tometric selection process. The w'BVRl filters are naturally suited to select Lyman break 
objects at z ~ 3 and are distinct from the UnGTZ filter set of Steidel & Hamilton (1992) in 
that they have the advantage of being readily available at most observatories. 

This marks the second step in an observing program designed to measure the DLA bias 
by searching for galaxies in emission that are clustered with DLAs or directly responsible for 
the damped absorption. Previous work by Gawiser et al. (2001) and Prochaska et al. (2002) 
detail the acquisition and analysis of galaxies associated with DLAs at z ~ 4 using deep BRI 
imaging. The ^ ~ 4 regime is more challenging due to a relatively low LBG number density 
to R^B < 25, a large number of \ow-z interlopers having similar colors, and the fact that the 
relevant interstellar absorption features are redshifted to wavelengths that coincide with a 
dense array of bright night-sky emission lines. By contrast, at z ~ 3 there are three times as 
many LBGs at a given magnitude, fewer low- 2; objects with similar colors, and the expected 
Lyman a and metal absorption line features are redshifted to a region with few night-sky 
emission lines and high CCD quantum efficiency. 

This paper is organized as follows: § 2 describes the deep imaging observations and 
reduction process, § 3 discusses the photometry and the photometric color selection, § 4 
describes the spectroscopic observations and analysis, § 5 addresses the use of spectroscopi- 
cally identified galaxies to refine the selection criteria, and § 6 provides a brief summary. The 
DLA-LBG cross-correlation and analysis are presented in Paper II. All magnitudes reported 
in this paper are given using the ABq^ (Fukugita et al. 1996) scale unless otherwise stated. 
We adopted the following transformations from the standard Johnson-Cousins magnitude 
system to the AB^^ system: u'j^^ = ^JC + 0.70, B^^ = Bjc - 0.15, Yab = ^ jc - 0.01, Rab 
= Rjc + 0.18, and Iab = he + 0.43. 

2. IMAGING OBSERVATIONS 

We chose QSO fields that contained at least one DLA in the redshift path consistent 
with the LBG detection volumes probed by ■u'BVRl photometric color selection. We also 
maintained a preference for systems near z — 3.0 where the photometric selection is most 
efficient and avoided systems with less than 3000 km sec~^ separation from the background 
QSO. Fields were chosen by their coordinate accessibility and, when more than one field was 
observable, we chose the field that would be observed through the lowest airmass and in the 
direction of lowest Galactic reddening. These considerations minimized the attenuation of the 
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light from the LBGs and eased their detection using reasonable integration times, especially 
in the w'-band. Prior to this survey, there were 10 known z ~ 3 DLAs in the selected 9 QSO 
sight lines. Our spectroscopy revealed a second DLA in the field PSS0808+5215 at z=2.936 
and was confirmed using the moderate resolution Keck ESI data of Prochaska et al. (2003). 
Table 1 lists the 9 QSO fields and 11 ^ ~ 3 DLAs used in this survey. Our m'BVRI color 
selection technique does not efficiently detect LBGs near the 4.18 DLA near PSS0957+3308. 
Although the imaging depths of the PSS0957+3308 field allowed for reasonable color selection 
of ^ ~ 4 LBG candidates using BVRI color criteria similar to those described in Gawiser et 
al. (2001), no follow-up spectroscopy has been acquired to date for these candidates. 

The imaging data were acquired during the interval 2000 April through 2003 November 
at the Palomar and Keck Observatories. We used the Carnegie Observatories Spectrograph 
and Multi-Object Imaging Camera {COSMIC; Kells et al. 1998) on the 200" Hale tele- 
scope at Palomar to obtain deep images for the first two fields in our survey. Both the 
PSS1432+3940 and JVAS2344+3433 fields were imaged in [/„BR filters using COSMIC 
equipped with a Tektronix 2048x2048 thinned CCD (0."286 pixel"^). Subsequent observa- 
tions in V and I were obtained using the Low Resolution Imager and Spectrometer {LRIS; 
Oke et al. 1995). The COSMIC field of view (FOV) is ~ 9.7 x 9.7 arcmin whereas LRIS 
FOV is ~ 5x7 arcmin. Two pointings with LRIS provided ~ 70% coverage of each 
COSMIC field. Table 2 presents a journal of the observations for the 9 fields in the survey. 

All other fields in our survey were imaged wholly using LRIS with u' , (f/„),B,V,R,(Rs), 
and I filters. LRIS is a dual arm instrument with both red and blue sensitive sides. Prior 
to January, 2001, all images were taken using the red arm equipped with a back-side illu- 
minated Tektronix 2048x2048 thinned CCD (0."215 pixel"^). Prom January, 2001 through 
June, 2002, C/„ and B-band images were obtained using the blue arm of LRIS with a SITe 
2048x2048 thinned CCD (0."211 pixel"^). After June, 2002, u' and B-band data were 
obtained using LRIS equipped with 2 Marconi 2048x4096 blue-sensitive CCDs (0."135 
pixel"^). 

Figure 1 plots the transmissions of the filters used in this survey. The transmissions 
of the Un and R^ filters are indicated by the dashed curves. The u' and C/„ filters were 
both used during the course of this survey as a result of their availability at Keck. Since 
the u' and Un filters have nearly identical bandpasses, effective wavelengths, transmission 
efficiencies, and application in this paper, the Un filter will be referred hereafter as the u' 
filter for simplicity, unless otherwise noted. The R^ filter was used wholly or partially to 
image three fields as the survey progressed in an attempt to increase the number of Lyman 
break candidate detections since it has a higher throughput and sharper transmission cutoffs 
than the Rjc- 
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AU observations were taken as close to the meridian as possible in an effort to minimize 
the effects of atmospheric extinction and dispersion. We imposed a ~ 20" non-repeating 
dithering pattern between consecutive images to minimize the effects of bad pixels and 
to allow for the creation of a sky flat by combining unregistered images (super-sky flat). 
Dithering by this amount had the effect of reducing the field size common to all images by 
~ 1' in both dimensions. The integration times were varied depending on filter, observing 
conditions, and telescope/instrument configuration in an attempt to achieve our targeted la 
magnitude sky fluctuation depths of it^^ = 29.0, B^^, ^ab, and Rab — 28.5, and Iab — 
27.8 with a seeing element goal of < 1". The average observed la magnitude sky fluctuation 
depths per seeing element were it'^^ = 29.1, B^^ = 28.7, Yab — 28.6, Rab — 28.7, and Iab 
— 27.6 with an average seeing element of 0."97 ± 0.2. Since the seeing for most fields was 
consistent, any u', B, V, or 1 image with a seeing FWHM differing from the corresponding 
R filter image by more than 0."2 was not used. This resulted in compatible seeing for each 
field in the final stacked images. For the three fields LBQS0056+0125, PKS0336-017, and 
PC1643+4631A, the stacked images consisted of only four filters. This was due to these 
fields having either a broad seeing FWHM from poor seeing or non-science grade images in 
a particular filter. These fields were only used in specific photometric selection and LBG 
candidate density results. 

The data were reduced using standard IRAF and IDL routines. Images were fiat- 
fielded using super-sky fiats whenever possible. In a few cases, the instrument condition or 
the particular night's agenda forced the use of a twihght fiat or dome fiat. Final images were 
combined with an in-house package that inversely weights each image by the measured vari- 
ance, scales by the exposure time, and creates la error images based on Poisson statistics. 
Photometric calibrations were calculated for each night by observing multiple photometric 
standard stars (Landolt 1992) in each filter over a full range of airmass. Conventional IRAF 
tasks, such as PHOT and F IT P ARAM S were employed to solve zero-point magnitudes, 
airmass correction coefficients, and appropriate color correction coefficients for the photo- 
metric data. Typical rms magnitude errors were ~ 0.04, ~ 0.02, and ~ 0.01 respectively. In 
the few cases where a unique solution was not found, we adopted fixed values for the air- 
mass correction coefficients. These values are: 0.56, 0.35, 0.17, 0.13, 0.07 for li'BVRI filters 
respectively at Palomar and 0.41, 0.19, 0.12, 0.11, and 0.07 for m'BVRI filters respectively 
at Keck. Finally, the photometric zero-points were corrected for Galactic extinction using 
E{B — V) values measured from the far-lR dust emission maps of Schlegel et al. (1998). 
The E{B — V) values for each field are given in Table 2. To correct our fields for Galactic 
extinction, we adopted the following relations interpreted from Cardelli et al. (1989): 

A{u') = 4.8 E{B - V), A{B) = 4.1 E{B - V), A{V) = 3.1 E{B - V) 

A{R) = 2.3 E{B - V), A{I) = 1.5 E{B - V) 
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3. PHOTOMETRY 

We defined our source detection catalogs in the R-band (Rg in 2 fields) after consid- 
ering the spectral profile of 2; ~ 3 LBGs and the imaging sensitivity of the LRIS and 
COSMIC cameras. An additional consideration was that ~ 30% of the PSS1432+3940 
and JVAS2344+3433 fields (imaged with the larger FOV of the COSMIC camera) had 
only ■u'BR filter coverage. The sources in all nine QSO fields were identified with the soft- 
ware package SExtractor (v. 2. 1.6) created by Bertin & Arnouts (1996). We varied the 
'(7 detection threshold' and 'minimum pixel area' input parameters by comparing the 
output segmentation map of each stacked R image with that of the segmentation map of 
the negative stacked R image to minimize spurious source detections. Typical SExtractor 
values were 1.6 cr detection threshold and 6 pixel minimum area. An in- house code was used 
to calculate the magnitudes and errors of each object detection in each filter. This code 
determines the sky flux by taking the bi-weight (Beers et al. 1990) of the sky in an annulus 
around each object containing at least 10^ pixels not flagged as an object in the segmentation 
map created by SExtractor. Magnitude errors were derived from variance in the sky flux 
and the systematic uncertainty in the photometric calibration. 

A goal in imaging the survey flelds was to achieve photometric completeness to the 
LRIS spectroscopic hmiting magnitude of R~ 25.5. The number of objects detected in our 
survey versus their R magnitude is shown in Figure 2. The upper histogram displays the 
number counts of sources for each filter in 0.2 magnitude bins for the five fields with it'BVRI 
imaging. A logarithmic fit to the binned number counts to R= 25.5 is used to estimate 
the completeness of our sample and is shown in the central plot. The lower plot displays 
the photometric completeness for all 5 filters. These fields are complete to R~25 and are 
> 80% complete to R=25.5 using this form of analysis. From the deep galaxy survey of 
Steidel et al. (2003), the average number of faint galaxies in the apparent magnitude range 
of 22.5 < TZ < 25.0 was found to be 25.9 with a field-to-field scatter of ±1.9 arcmin"^. 
We report an average faint galaxy number of 24.9 ± 3.5 arcminute"^ that meet our source 
identification criteria over the same magnitude range in the R-band. 

To facihtate the design and milling of reliable spectroscopic multisht masks, accurate 
astrometry was performed over the full area of each field using USNO A2.0 and USNO Bl.O 
(Monet et al. 1998, 2003) astrometric catalogs. We then used the PAL_COORDS routine^ to 
determine the astrometric solution for both LRIS and COSMIC cameras. Typical position 
errors of ~ 0."2 were reported by this method and confirmed when fine aligning the slitmasks 
on sky. 



^Available by anonymous ftp at: ftp.astro.caltech.edu/palomar/paLcoordinates/paLcoordinates.doc 
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3.1. Photometric Color Selection 

We chose to identify the color of all LBG candidates in a consistent manner from their 
coarse spectral features. We searched for star forming galaxies having a strong continuum 
long- ward of rest-frame Lyman a (1215A) and a flux decrement short-ward of the rest-frame 
Lyman limit (9 12 A). This decrement is due to absorption by H I optically thick at the Lyman 
limit which is intrinsic to the galaxy and present in foreground systems. At z ~ 3, these 
spectral characteristics are redshifted to optical wavelengths. The photometric selection of 
2; ~ 3 LBG candidates in color-color space has been spectroscopically verified with ground- 
based telescopes (Steidel et al. 1996a) using UnGlZ broad band filters and with the Hubble 
Space Telescope [Steidel et al. (1996b) and Lowenthal et al. (1997)] using similar bandpasses. 
Here we use five widely available filters, u\ B, V, R, and I, to formulate our LBG candidate 
selection criteria. 

Galaxies meeting our selection criteria exhibit virtually no fiux in the tx'-bandpass, a 
varying flux decrement in the B-bandpass due to the Lyman a forest, and relatively constant 
flux in the VRI bandpasses. Figure 3 displays a graphical representation of this in the u' , 
B, and V fllter bandpasses. The placements of the Lyman a and Lyman limit breaks in 
the continuum are readily apparent. We derived theoretical galaxy colors using six star- 
forming galaxy templates (Yahata et al. 2000) and simulated the color-color evolution of 
these galaxies to high redshift. Figure 4 plots the expected {u' — B) vs. {B — R), [u' — V) 
vs. (F — i?), and {u' — V) vs. {V — /) colors of a starbursting galaxy (dashed green curve) 
and an elliptical galaxy (dotted blue curve) as they are evolved from z = Q. Star symbols are 
used to denote the colors of Galactic stars from the Gunn & Stryker (1983) catalog. From 
approximately 2.5 < z < 3.5 in each plot, the starbursting galaxy is indicated by two solid 
red lines. The separate red lines indicate the decrements in the continuum of the galaxy 
caused by a reasonable range of Lyman a forest absorption (Madau 1995). The expected 
colors of these objects were used to determine the appropriate color criteria that would 
efficiently detect Lyman break candidates at 2; ~ 3 and reject low redshift objects. The color 
selection regions are shown as the boxed regions in the upper-central and upper-left portions 
of the plots. 

Objects were selected as LBG candidates that satisfied the following constraints: 

{u' - B)ab > 1.1, (4) 

{u' -V)ab>I.Q, (5) 
0.6 < (S-i?)^B <2.1, (6) 
{V - R)ab < 0.6, (7) 
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{V - I)ab < 0.6, 



(8) 
(9) 



20.0 < Rab < 25.5. 



LBG candidates for the first two spectroscopic slitmasks of the field JVAS2344+3433 
were determined using the {u' — B) vs. {B — R) color-color plot only (constraints 4, 6, & 9). 
This was necessary prior to complementary V and I imaging at Keck and remains the color 
selection criteria for the regions of fields PSS1432+3940 and JVAS2344+3433 having only 
m'BR imaging. LBG candidates in fields where the images in all five filters were usable 
were required to meet the full color criteria (constraints 4 — 9). Constraints 4, 5, and 6 
were imposed to reflect the relative discontinuities in the continuum at rest-frame 912A and 
12 15 A from the continuum at rest- frame ~1300A. Constraints 7 & 8 were used to ehminate a 
large number of low redshift reddened eUiptical galaxies and Galactic stars. The R magnitude 
(equation 9) of all objects was constrained on the bright end (Rab > 20.0) to avoid saturated 
sources in the combined images and constrained on the faint end {Rab < 25.5) to only include 
objects bright enough to meet the spectroscopic limit of LRIS on Keck. The imposition of 
color constraints similar to the ones mentioned above selects the majority of the galaxies at 
z 3 from field objects but does not allow for the selection of 2; ~ 3 objects with excessive 
intrinsic far-UV extinction. Therefore, this survey is not complete to all possible types of 
galaxies at 2; ~ 3, but does include a significant number of galaxies to provide meaningful 
statistics toward the goal of cross-correlating LBGs with DLAs. 

Objects meeting the full m'BVRI color criteria (constraints 4 — 9) were given the highest 
priority as slitmask targets. Objects meeting only u'VRI and u'BRl color criteria were given 
successively lower priorities. To compensate for photometric errors due to both statistical 
and systematic uncertainties, objects within 0.2 magnitudes outside of these constraints were 
also identified as LBG candidates and given respectively lower priorities. LBG candidates 
within a separation of ~ 10" in radius from a QSO sight-line were marked and chosen over 
other objects that were in conflict in the dispersion direction of the particular slitmask in an 
effort to identify the galaxy responsible for the damped absorption. A separation of ~ 10" 
represents the maximum impact parameter of an absorber inferred from absorption line 
statistics (Storrie-Lombardi & Wolfe 2000). In some cases, the conflicting LBG candidates 
were placed on subsequent shtmasks. No other factor contributed to the choice of candi- 
dates for spectroscopy except the physical constraints of the shtmask itself, as described in 
§ 4. The three color-color plots in Figure 5 display the data from 5 flelds (PSS0808-I-5215, 
PSS0957+3308, BRI1013+0035, PSS1057+4555, and PSS1439+3940) having full m'BVRI 
color selection. In this figure, all 20.0 < R < 25.5 objects in the five fields are shown as 
red points. Green diamonds are overlayed to indicate those objects that met the imposed 
m'BVRI color selection criteria (including objects within 0.2 magnitudes) and the black as- 
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terisks denote the objects that have secure z >2 spectroscopic redshift identifications. 

Overall, a total of 796 objects were photometrically selected using these color criteria. 
In the four fields where all five filters could not be utilized, candidates were chosen that met 
the color criteria of the usable filters. Firstly, the initial shtmasks for the JVAS2344+3433 
field included LBG candidates chosen by ■u'BR color selection prior to V and I filter imaging. 
Secondly, the LBG candidates in both the LBQS0056+0125 and PC1643+4631A fields were 
chosen using -u'VRI color selection. An instrument failure resulted in non-science grade B 
filter images for the LBQS0056+0125 field and time constraints forced the full bandpass 
of the B filter to be compromised by a dichroic centered at 4600A for the PC1643+4631A 
field. Lastly, the LBG candidates in the PKS0336-017 field were chosen using w'BRI color 
selection due to comphcations from merging a broad and shallow FWHM of l."7 in the V 
filter. All LBG candidates in this field were required to meet an additional color criterion 
of: 

0.6 < < 2.1. (10) 

Table 3 displays the area common to all filters for each targeted field and summarizes 
the number of LBG candidates detected using w'BR(I), m'VRI, and m'BVRI color selection 
criteria and the resulting m'VRI and m'BVRI LBG number densities. In the five fields with 
m'BVRI photometry, the two-dimensional LBG color candidate density is 0.79 arcmin~^±0.25 
arcmin"^. Comparison of the average LBG candidate density with that of similar surveys 
can be used as a photometric selection check. Prom Steidel et al. (2003), the LBG candidate 
density using UnGTZ I filters is 1.6 ±0.5 arcmin"^ from 2347 selected objects in 17 fields. Our 
average LBG candidate density using the relatively comparable u'VKl filters is consistent 
at 1.45 ±0.66 arcmin"^ from 365 objects in 6 applicable smaller fields. Although the zx'VRI 
criteria select a larger number of objects than the n'BVRI criteria, we find the w'VRI method 
to be less efficient in resulting z > 2 spectroscopic confirmations (see § 5). 



4. SPECTROSCOPIC DATA 

To optimize spectroscopic slitmask design, objects on the LBG candidate lists were 
visually inspected in the two-dimensional images to eliminate false detections (e.g. satellite 
trails, diffraction spikes, asteroid tracks, and CCD defects). LBGs have a half-light radius 
of < 0."5 from Hubble Space Telescope images (Steidel et al. 1996b; Lowenthal et al. 1997; 
Gardner et al. 2001) and should appear nearly point-like with the resolution of LRIS and 
COSMIC. Therefore, several objects were rejected if they were overly extended [although 
some z ~ 2 objects appear to be extended (i.e. Erb et al. 2004)]. Each slitmask was then 
designed to include the highest priority objects and to have empty regions filled with objects 



- 12 - 



of successively lower priorities. We chose to use 6" minimum slitlet lengths to maximize the 
number of objects observed without compromising sky subtraction. Slitmasks were designed 
using the software AUTOSLIT 3.07 Typically, each slitmask was able to include 25 — 35 
objects which resulted in an average slitlet length of ~ 16" for each candidate. 

Multislit spectroscopic observations of the LBG candidates were obtained between 2000 
November through 2004 February, using both the red arm and new blue arm of LRIS on the 
Keck 1 telescope. We used the 300 line mm~^ grating (red side) and the 300 line mm~^ grism 
(blue side), both blazed at 5000A in first order. The red side grating has a dispersion of 
2.47A pixel^^ over the 2048 x 2048 CCD and was tilted to allow a minimum wavelength range 
of 3800 — 7OOOA for each slitlet regardless of their placement on the mask. The blue side grism 
has a dispersion of 1.43A pixel"^ (15yum pixel size) over two Marconi high quantum efficiency 
(QE) blue sensitive 2048 x 4096 CCDs. The grism is fixed and designed to allow a possible 
instrumental wavelength range of 2210 — 8O6OA for multislit spectroscopy. This results in 
an observational minimum wavelength coverage from the atmosphere cut-off at ~ 3200 to 
~7500A for each slitlet. We designed the widths of the slitlets on most slitmasks to be 
l."5 to account for errors in the astronietry, atmospheric dispersion over long integrations, 
and moderate seeing. Since these are sky-dominated objects, when possible, we designed 
alternate slitmasks with narrower slitlet widths (l."0 — l."3) to accommodate good seeing 
and observations at low airmass through the meridian. The spectral resolution evaluated 
from night sky emission lines was ~12A FWHM and, because Lyman break objects are 
near point sources, typical object resolution was seeing limited and resulted in an effective 
spectral resolution of ~9A FWHM. 

The total exposure time for each sfitmask was 7200 seconds when using the red arm and 
~ 3600 — 4800 seconds when using the blue arm of LRIS. The higher throughput, lower 
read-out noise, and higher CCD quantum efficiency of the blue arm permitted considerably 
shorter exposure times to achieve equivalent magnitude depths as compared to the red arm. 
Each slitmask was observed in a series of consecutive 1200 — 1800 second exposures with 
small dithers (0."5 — l."0) along the direction of the slits to sample different parts of the 
CCD, minimize fringing effects, avoid hot pixels, and provide better sky subtraction options. 
The details of the spectroscopic observations are summarized in Table 4. The slitmasks were 
designed to be observed at position angles equaling the parallactic angle determined mid- 
way through their scheduled time of observation. In some instances (to help compensate for 
uncertain weather and instrument conditions), two masks were designed with position angles 
reflecting two different times of possible observation. This practice minimized wavelength 



^Autoslit and support files are available at littp://astro. caltech.edu/~pls/autoslit/ and 
ftp :// astro . caltech . edu/users/pls/ autoslit / 
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dependent flux loss through the slit due to atmospheric dispersion as the slitmasks were 
only observed within 30 minutes of their scheduled time. Calibration exposures were taken 
either before, mid-way, or after the science exposures to minimize the effects of instrument 
flexure. Flat field exposures were obtained using an internal halogen lamp. Internal arc lamp 
exposures used for wavelength solutions were taken with a combination of Cd, Zn, Hg, Ne, 
Ar, Kr, and Xe, depending whether the science exposures were taken with the red or blue 
arm of LRIS. For proper flux calibration, two or more spectroscopic standard stars (Massey 
et al. 1988) were observed each night at different airmass to rephcate the range of airmass 
sampled by the slitmasks during the course of the night. 

We used IRAF and in-house IDL scripts to perform the initial reduction of the CCD 
data including overscan correction, gain normalization, trimming, bias subtraction of all 
images. We used the IRAF routine bogus. cl ^ to handle the division of each full two- 
dimensional multislit image into separate two-dimensional slitlet images of each object and 
for the application of calibration frames. The bogus . cl routine calls IRAF tasks for flat 
fielding the individual slitlet images by their corresponding continuum-normalized portion 
of the fiat field image. Background sky subtractions were determined by using low-order 
polynomial fits appropriate to the length of each slitlet to background sky sections. The 
routine also performs fringe subtraction, identifies cosmic rays morphologically, and shifts 
and stacks the spectra. Additionally, bogus . cl creates two-dimensional shtlet images during 
selected steps for later inspection and analysis. We chose to operate interactively on all 
spectra to achieve our goals considering the faint magnitude of typical objects in our survey. 

We traced the reduced two-dimensional stacked spectra using IRAF spectroscopic anal- 
ysis routines and fitted low-order (usually first-order) polynomials to subtract any remain- 
ing background sky offset. Corresponding arc lamp traces of each object were extracted in 
the identical manner as the object. Wavelength solutions were obtained on the extracted 
arc lamp spectra using low order polynomials with residuals < lA. Further refinement to 
the wavelength solution was performed by comparing the arc lamp wavelength solution to 
multiple night sky emission lines corresponding to each object and applying the calculated 
sub-pixel offset. Each spectrum was then corrected for extinction, fiux calibrated using sen- 
sitivity functions obtained from the corresponding spectrophotometric standard stars, and 
corrected to vacuum wavelengths. 



bogus. cl is available at: littp://zwol£k;inder.jpl. nasa.gov/~stern/homepage/bogus. html 
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4.1. Spectroscopic Identification and Analysis 

LBGs display identifiable rest-frame UV spectral features from 912A to ~1700A. For 
galaxies at 2; ~ 3, these features are shifted to optical wavelengths from ~ 3200 — 7500A. 
The objects in this survey were identified mainly using Lyman a and interstellar emission 
and absorption features described by Pettini et al. (2000), Adelberger et al. (2003), and 
Shapley et al. (2003). These features have typical equivalent widths of ~ 5 — lOA and can 
be relatively easy to identify even though the spectra have a signal-to-noise ratio of only a 
few. 

We chose to separate the LBGs detected by this survey into two simple populations: 
(a) emission-identified LBGs - galaxies showing Lyman a in emission and (b) absorption- 
identified LBGs - galaxies showing Lyman a in absorption. Confirmation for all LBGs 
involved the individual search through each continuum for a possible decrement due to the 
Lyman a forest short-ward of the Lyman a feature at rest-frame 1215A, a second decrement 
(where spectral coverage allowed) due to Lyman limit absorption short-ward of rest-frame 
912A , and the identification of expected interstellar absorption and/or emission features. 
The differing coarse spectral characteristics of LBGs as described by Shapley et al. (2003) 
prompted slightly different identification expectations for each population. Spectra exhibit- 
ing Lyman a in emission and appropriate continuum breaks were additionally inspected for 
the presence of interstellar absorption lines (usually weak) and a fiat continua long-ward of 
Lyman a to aide in their identification. Spectra with Lyman a in absorption were inspected 
for the appropriate continuum break characteristics and multiple strong absorption lines 
with a somewhat reddened continua long-ward of Lyman a. At no time were the spectra 
required to have all of these features, although nearly all identified galaxies fulfilled these 
expectations. 

In practice, and depending on the wavelength coverage of an individual spectrum, usu- 
ally 6 or more absorption lines were used to secure the redshift of a particular LBG. Due 
to the overall low signal-to-noise nature of these spectra, it becomes necessary to assign 
secure and probable LBG redshifts. Both the secure and the probable LBGs exhibit the 
characteristic continuum profile. Relevant expected stellar and interstellar metal fines used 
for identification included: Si II 1193, 1260, 1304, 1527, Si IV 1394, 1403, C II 1335, C IV 
1548, 1551, O I 1302, Fe II 1608, He II 1640, Al II 1671, Al III 1855, 1863 , and GUI] 1909 
among others. Low redshift galaxies were identified by comparison to templates of E-I-SO, 
E-l-A, and starbursting galaxies. Since the break in the continua of low redshift galaxies 
near 4OOOA can mimic the Lyman break of LBGs, low signal to noise absorption spectra 
must be carefully checked to avoid misidentification. Emission lines, such as [O II] 3728, [O 
III] 4959, 5007, [N II] 6584 and [S II] 6724 and typical absorption and/or emission features 
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such as the Hydrogen series, Ca II H & K hnes, Mg II 2796, 2804, Mg I 2852, 5175, He I 
4471, 5876, 6678, 7065 and Na I 5890, 5896 were used to confirm low redshift galaxies and 
Galactic stars. 

In total, we acquired spectra of 529 objects from the 796 LEG candidates that met either 
u'BR(V), w'VRI, or it'BVRI photometric selection criteria. Using an in- house analysis routine, 
we interactively examined each one-dimensional and two-dimensional spectrum individually 
for identifiable features following the prescription described above. All objects were then 
placed in one of the following seven categories: (0) stars, (1) secure z < 2 galaxies, (2) 
probable z < 2 objects, (3) probable z > 2 objects, (4) secure z > 2 galaxies, (5) high 
redshift QSOs, and (6) unidentified objects. Objects placed in categories 2 and 3 were 
identified as low or high redshift objects respectively, but did not satisfy the identification 
requirements as rigidly as categories 1 and 4. The probable identifications usually resulted 
from lower signal-to-noise and lack of apparent emission lines yet still showed reasonable 
continuum breaks and the tentative identification of several to many absorption lines. Some 
objects were placed in the probable LEG category due to a slight ambiguity in redshift. 
Objects were classified as unidentified (category 6) if they displayed very low signal-to-noise 
(possibly from a portion of these objects falling outside the slitlet) or a strong ambiguity in 
redshift assignment. These are listed for each field in Table 5 with the separate categories 
indicated by the parentheses. 

Two composite LEG spectra are shown in Figure 6. Each spectrum was constructed 
by stacking the rest-frame spectra of ~ 20 individual LEGs that exhibit similar properties. 
The upper spectrum is a composite of LEG spectra with prominent Lyman a emission. This 
spectrum displays visibly weak (narrow) interstellar metal absorption lines and a relatively 
fiat continuum long- ward of 12 15 A. The Lyman a feature is redshifted with respect to the 
interstellar hnes by an average of ~ 600 km sec~^, in agreement with Adelberger et al. 
(2003). This is interpreted to be due velocity differences within the galaxy arising from 
stellar and sup ernovae- driven winds. The higher redshift of the Lyman a feature results 
from dust absorption of the resonantly scattered photons in the blueshifted portion of the 
expanding galactic shell. Emission-identified LEGs comprise ~ 46% of our sample. This is 
a greater fraction of LEGs of an equivalent type than found by Steidel et al. (2003) and may 
be due to a more conservative redshift identification process. The lower spectrum combines 
LEG spectra having Lyman a in absorption, revealing a somewhat reddened continuum and 
noticeably stronger (broader) interstellar absorption lines. The assessment of the overall 
continuum profiles of LEGs from our data is in agreement with the results of the analysis 
of intrinsic LEG properties by Shapley et al. (2001, 2003). The fiuxing of the spectra long- 
ward of ~1500A in Figure 7 is not reliable due to the fact that no order blocking filter was 
used for these observations. Since standard stars used for fiuxing tend to be very blue, their 



-16- 



second order light can significantly affect the spectrum at observed wavelengths long-ward 
of ~6400A. This corresponds to rest-frame wavelengths long-ward of ~1600A for objects at 
z 3. Moreover, the composite spectra comprise LBGs of varying redshifts. As a result, the 
effects from the flux degradation is spread over rest-frame wavelengths greater than ~1500A. 

To test our identification process for any bias toward magnitude, we performed a two- 
sided Kolmogorov-Smirnov (K-S) test of the secure and probable LBG R magnitude distri- 
butions. This resulted in an 88% K-S confidence level for the null hypothesis, making it very 
unlikely that they are drawn from separate distributions. Figure 7 shows a histogram of the 
secure and probable LBGs versus R magnitude binned in AR=0.2. Additionally, we tested 
for any bias in magnitude toward LBGs displaying Lyman a in emission versus LBGs dis- 
playing Lyman a in absorption since the presence of a strong emission feature could enhance 
the confidence of identification. This is especially the case for low signal-to-noise spectra. 
Figure 8 presents R magnitude histograms of emission-identified LBGs with absorption- 
identified LBGs. The K-S tests on the secure, probable, and combined LBG distributions 
resulted in a 29%, 46%, and 51% confidence for the null hypothesis, respectively. In all three 
cases, the values are consistent with the hypothesis that the secure emission-identified and 
absorption-identified LBGs are drawn from the same parent distributions. 

The ability to securely identify an emission-identified LBG relative to absorption-identified 
LBG for the faintest objects tends to diminish with magnitude but will actually increase 
beyond R~ 25. This is reflected in Figure 8 where it can be seen that there are no emission- 
idcntiflcd LBGs and 15 absorption-identifled LBGs with R > 25. Although this is partially 
due to the presence of the Lyman a feature adding confidence to these lowest signal-to-noise 
spectra, it is also an effect of the photometric selection criteria which is designed to select 
galaxies based on the decrements in their continua. Due to the finite magnitude depths of the 
u' and B-band images combined with strict {V — R) and {V — I) constraints, only the bluest 
objects selected with R ~ 25 and greater can meet the selection criteria. The bluest LBGs 
are typically emission-identified LBGs. Because of this, the 15 R> 25 absorption-identified 
probable LBGs are the least confident of the sample. 

The redshift distribution of the identified LBGs is another regime to test our identifica- 
tion process and lends insight into the overall distribution of LBGs detected by this selection 
method. Figure 9 presents histograms that display secure, probable, and a combination 
of the two LBG categories in the upper, central, and lower histograms respectively. The 
broad agreement in the absorption-identified LBG redshift distributions lends confidence 
to the overall redshift assignments. The separate sets of emission-identified and absorption- 
identified LBGs show a small offset in the centers of their distributions. The secure emission- 
identified LBG distribution centers at ^ = 3.07 ± 0.33 while the secure absorption-identified 
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LBG distribution centers ai z — 2.92 ± 0.25. Nevertheless, a K-S test on the two secure 
distributions resulted in a 21% confidence of the null hypothesis and is consistent with the 
two distributions being drawn from a similar parent distribution. 

Performing the K-S test on the probable LBG redshift distributions resulted in a 99% 
confidence that the emission-identified and absorption-identified LBGs are likely to be drawn 
from different distributions. This assessment is considered in future analysis and can be 
seen by the apparent lack of probable emission-identified LBGs with z < 2.9. Although 
emission features add to the general confidence of spectral identifications, emission-identified 
Lyman break objects at the higher end of the redshift range probed by our method become 
less confident since the majority of their expected stellar and interstellar features begin to 
compete with a dense array of night sky emission lines. This is also the case for the most 
common low redshift interlopers. The LRIS instrument is less capable of subtracting these 
lines effectively for very low signal-to-noise spectra. 

Prom the full spectroscopic sample, we arrive at a secure (category 4) LBG number den- 
sity of 0.37 arcmin"^ with a field-to-field scatter of ±0.21 arcmin"^ and 0.56 ±0.27 arcmin"^ 
for the combined set of secure and probable LBGs (categories 3 and 4). Approximately two- 
thirds of the photometrically selected R< 25.5 LBG candidates have spectroscopic coverage 
and approximately 63% of the obtained spectra are identified. Assuming similar results 
for the remaining R< 25.5 LBG color candidates gives a projected LBG number density 
of 0.78 ± 0.46 arcmin"^ for the secure LBGs and 1.21 ± 0.49 arcmin"^ for the secure and 
probable LBGs. 

4.2. The Photometric Selection Function 

The selection criteria described in § 3 were designed to be sensitive to objects over a 
redshift path from 2.6 < z < 3.4. The previous large LBG survey at z ~ 3 by Steidel et al. 
(2003) probes a similar redshift path with a distribution of (z) = 2.96 ±0.29 for their catalog 
of 940 identified spectra. They used their spectra to establish several fundamental properties 
of LBGs at 2; ~ 3. Our color selection technique was intentionally designed to select LBGs 
over the same redshift path to utilize the LBG properties defined by Steidel and coworkers. 
The redshift distributions of the LBGs detected by this survey are presented in the upper 
two histograms of Pigure 10. The top histogram displays the 137 secure LBG redshifts from 
2.4 < z < 3.6 binned in Az = 0.04. The central histogram represents the distribution 
of the combined set of LBGs, which includes 74 additional probable LBG redshifts. The 
secure LBG distribution has a mean of (z) = 3.00 with a standard deviation of ±0.32. The 
combined set of 211 objects is nearly identical with (z) = 3.02 ± 0.32. These distributions 
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include the five faint 2; ~ 3 QSOs discovered in this survey. 

Observed redshift distributions combine the actual LBG redshift distribution with the 
effects of the photometric selection method. To measure the clustering properties of LBGs 
associated with DLAs, the observed LBG densities per redshift interval must be corrected 
for the photometric selection effects caused by the efficiencies and transmissions of the in- 
struments and filters as well as the photometric uncertainties. These corrections can be used 
to recover the true LBG background redshift distribution from the observed LBG redshift 
distributions once the true photometric selection function is determined. Both observed 
LBG redshift distributions from the data (Figure 10) were fit by Gaussian functions (solid 
curves) using the mean, standard deviation, and the same total number of objects of each 
distribution. The normahzed residuals to the fits with respect to the expected Poisson fiuc- 
tuations per bin are —0.05 ± 0.95 (0.03 ± 0.89) for the secure (combined) set of LBGs and 
are over-plotted on Figure 10 as dashed curves. The fits indicate that in both cases the 
product of the photometric selection function and the true LBG redshift distribution are 
well described by Gaussian functions. 

To model the true photometric selection function, we generated a catalog of 1000 star- 
bursting galaxies with random magnitudes drawn uniformly from 23.5 < R < 25.5 and 
random redshifts from 2.4 < z < 3.6. These galaxies were then assigned respective u'BVl 
magnitudes according to their assigned redshifts using the theoretical galaxy templates 
of Yahata et al. (2000). The quantum efficiency of the CCDs, atmospheric attenuation, 
and transmission characteristics were factored into the response of each filter. We fit a 
function to the photometric errors of the data with respect to magnitude and assigned 
Gaussian random errors multiplied by this function to each galaxy accordingly. Under the 
assumption that the LBG spatial density varies slowly with redshift, we assumed a uni- 
form background LBG redshift distribution and corrected the expected observed LBG back- 
ground to include the fact that a standard candle diminishes by ~ 1.1 magnitudes from 
z — 2.4 — 3.6 {h — 0.72, Qm = 0.3, fl\ = 0.7). In doing this, we incorporated the fact that 
our survey is magnitude limited by only considering objects with a corrected magnitude of 
R< 25.5. 

The results of this simulation are shown in the lower histogram in Figure 10. This 
simulation shows that the product of a uniform LBG redshift distribution and the simulated 
true photometric selection function is well described by a Gaussian (solid line) and results 
in a distribution with nearly identical mean and standard deviation as the data. In an 
identical manner as the data, the normalized residuals of the Gaussian fit to the simulation 
are over-plotted as the dashed curve and have a mean of —0.11 ± 0.76. 

Since the simulated redshift distribution which incorporated the effects of our photo- 
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metric selection method and assumed a uniform background LBG redshift distribution is 
shown to be described well by a Gaussian function, and since the observed redshift distribu- 
tions can similarly be described by Gaussian functions, it is reasonable to assume a uniform 
background number density from 2.4 < z < 3.6. Knowledge of the LBG background number 
density per redshift interval is required for proper determination of the DLA-LBG clustering 
amplitude. Although the true background LBG number density from 2.4 < z < 3.6 may 
not be strictly uniform, any slow varying non-uniformity in the large-scale distribution of 
LBGs will not affect the DLA-LBG clustering amphtude significantly. Indeed, the spatial 
distribution of LBGs is clustered on small scales, but bin-by-bin averaging over the 9 fields 
and 11 z ~ 3 DLAs in this survey will in effect "wash out" an individual field's clustering 
signal. 



5. COLOR SELECTION EFFICIENCIES AND CRITERIA REFINEMENT 

The u'BVRI photometric color selection technique presented here is an efficient means 
to detect LBGs over the intended redshift path of 2.6 < z < 3.4. Although the original 
intent of this survey was to perform it'BVRI photometric color selection in all 9 fields, 
various instrument and imaging concerns, described above, required the use of a subset 
of the fuH color criteria in four fields: LBQS0056+0125, PKS0336-017, PC1643+4631A, 
and JVAS2344+3433. Restricting these fields to a subset of the w'BVRl filters did not 
significantly hinder our ability to detect LBGs. For the five fields with m'BVRI imaging, 
the number density of LBG candidates meeting m'BVRI color criteria allowed room on each 
slitmask to include objects meeting lesser color constraints. This allowed an additional 
opportunity to examine the efficiencies of these alternate criteria. 

Table 6 summarizes the efficiencies of the m'BVRI color selection criteria and subset 
criteria from the spectroscopic data. Adhering to the imposed ■u'BVRI color selection cri- 
teria outhned in § 3.1 and the spectroscopic identification process in § 4.1, we find 79% of 
the identified objects selected for spectroscopy as secure and probable z > 2 Lyman break 
objects, 0% were Galactic stars, and 21% were secure or probable z < 2 galaxies. The 
efficiency is clearly influenced by the conservative nature of our redshift assignments and 
subjective object placement in the probable categories. Re-categorizing the probable spec- 
troscopic placements would lead to as much as 95% and as few as 59% of the identified 
objects with z > 2. 

Nearly 60% of the unidentified spectra can be attributed to cirrus during the 2003 
February observing run for the PKS0336-017 field and an instrument failure for the 2002 
December observations of the PSS1057+4555 and JVAS2344+3433 fields in which the CCD 
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shutter could not be closed. In the latter case, we corrected for the additional light striking 
the CCD during read-out by subtracting a one-second image taken immediately after each 
observation image. This was surprisingly effective in removing global flux levels, but was by 
no means perfect, especially in the complete removal of bright galaxy and star trails. 

A review of the placement of spectroscopically confirmed LBGs on the color-color 

plots allows an opportunity to assess and possibly refine the color selection criteria. Fig- 
ure 11 displays the spectroscopic results from the five fields (PSS0808+5215, PSS0957+3308, 
BR11013+0035, PSS1057+4555, and PSS1432+3940) that have m'BVRI color selection. 
Zoomcd-in views of the selection regions arc shown and indicated by the boxed regions 
(dashed lines). LBGs (both secure and probable) selected by their m'BVRI color criteria are 
shown as black diamonds, those selected by only their m'BR(I) color criteria are shown as 
blue squares, and those selected by only their u'YRl colors are shown as green triangles. All 
objects with z < 2 identifications are shown in red. The z < 2 galaxies and Galactic stars 
meeting ;i'BVRI color criteria are shown collectively as asterisks and those meeting m'BR(I) 
or u'VRI color criteria are shown as plus signs or 'X's respectively. 

The practice of obtaining spectra of the lower priority candidates revealed and helped 
quantify the fraction of the galactic population at 2; ~ 3 that is missed when using the 
full set of m'BVRI filters. Additionally, objects detected and confirmed using only t('BR(I) 
color criteria are missed when using only -u'VRl color selection and vice-versa. This can be 
seen in Figure 11. The bhic squares to the right of the selection region in the {u' — V) vs. 
{V — R) plot were detected and confirmed using w'BR(l) color selection. Similarly, the green 
triangles below and to the left of the selection region in the {u' — B) vs. (B — R) were detected 
and confirmed using u'VRl color selection. The original m'BVRI color criteria (dashed fines) 
provide a very efficient method to discern the largest possible number of Lyman break objects 
from fields objects with a minimum number of undesired low-redshift objects. To acquire a 
more complete sample of all populations of galaxies at ^ ~ 3, the data suggest refinements 
to the selection regions if either all or a subset of the u'BVRI filters are used. 

Prom the {u' — B) vs. {B — R) plot in Figure 11, it appears that the original {u' — B) 
color cut is effective and any relaxation of the criterion beyond 0.2 magnitudes would include 
a large number of less desired field objects (review Figure 5 for the field object distributions). 
Additionally, the extension of the {B — R) color cut to low values would help to include or 
confirm objects detected via -u'VRl color selection. Due to the lack of field objects in this 
region, this extension does not risk a loss to efficiency. A refinement to this plot (solid lines) 
would be to set the {u' — B) criterion to {u' — B) > 0.9 and modify the original color criterion 
from 0.6 < (B - R) < 2.1 to (B - R) < 2.4. 

A modification of the color criteria in the {u' — V) vs. {V — R) color-color plot (solid 
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lines) is intended to include objects that would be detected via m'BR(I) color selection. 
This refinement would modify the current criteria from {u' — V)> 1.6, {V — R) < 0.6 to 
{u' -V)> 1.2, {V -R)< 1.2, {u' -V)>IA{V - R)+ 0.5. Unfortunately, the distribution 
of low redshift interlopers extends into the modified selection region. This refinement would 
include a broader sample of galaxies at the risk of diminishing the LBG detection efficiency. 

Several objects chosen by only their m'BR(I) colors (blue squares), show iip to the right 
of the original selection region in the (V — R) and (V — I) color-color plot. To include 
these objects and to search for others with reasonable SEDs without the high cost of many 
interfering field objects, a suggested refinement to the {u' — V) vs. {V — I) color criteria 
is shown (solid lines). The suggested refinement involves modifying the original color cuts 
from {u' -V)> 1.6, {V - I) < 0.6 to {u' -V)> 1.2 and {u' -V)> 3.8 {V - I) - 2.6. 

Finally, in our analysis of the 529 color selected objects, we also analyzed as many field 
objects as possible that randomly fell into our slitlets. These objects provide a small sample 
of the full population of galaxies at high redshift free from color selection bias. Prom the 132 
serendipitous field objects, and using the categories described above, we found Galactic 
stars, 16 secure z <2 galaxies, 27 probable z <2 objects, 8 probable z >2 galaxies, 4 secure 
z > 2 galaxies, QSOs, and 77 unidentified objects. Three of the secure z > 2 galaxies are 
z > 3.65 and would not have been detected by our selection criteria. The fourth galaxy is an 
emission-identified LBG at 2; = 3.263 that falls below the spectroscopic magnitude cut-off 
for our survey of R=25.5. Of the probable z > 2 galaxies, 3 have z > 3.55 and the remaining 
5 were found to have 2.94 < z < 3.29 and R> 25.5. Although the 55 serendipitous objects 
identified sample a very small area of the sky relative to the entire survey, none of these 
objects were found within the redshift path probed by our selection method and brighter 
than our magnitude hmit of R=25.5. These statistics are consistent with the idea that our 
selection method detects the majority of objects at 2; ~ 3. 



6. SUMMARY 

This paper presents the first spectroscopic survey designed to detect galaxies associated 
with known DLAs at 2; ~ 3. A primary goal of this survey is to measure the spatial clustering 
of LBGs and DLAs to determine the DLA-LBG cross-correlation and infer the typical mass 
of DLAs in the context of CDM cosmology. Since DLAs may have a range of mass, a survey 
for LBGs in multiple QSO fields with known DLAs is required to achieve proper statistics. 
We present and analyze the imaging and spectroscopic data for objects in 9 QSO fields 
having 11 known 2; ~ 3 DLAs, including one DLA discovered by this survey, for the purpose 
of measuring the DLA-LBG cross-correlation. 
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We selected a total of 796 objects as LBG candidates using the color criteria described 
in §3 from 17,343 sources detected over the 465 arcmin^ in our survey. We find an LBG color 
candidate number density of 0.79±0.25 arcmin"^ in the five fields with it'BVRI imaging and 
1.45 ±0.66 for six fields with w'VRI imaging. 

We obtained spectroscopic observations of 529 LBG candidates and identified 339 red- 
shifts in which 211 are secure and probable z > 2 LBGs with a distribution of {z) = 
3.02 ± 0.32. Of the 211 LBGs, 46% display Lyman a in emission and the remaining display 
Lyman a in absorption. We find no bias in the confidence of our LBG identifications with 
magnitude and no bias in magnitude based on individual LBG spectral features. Addition- 
ally, we find no statistically significant difference in the redshift distributions between LBGs 
identified by spectral features yet we find an effect in confidence of our probable identifi- 
cations with redshift. We obtained spectroscopic data for ~ 67% of the photometrically 
selected LBG color candidates in nine fields and identified 63% of the obtained spectra. 
We find an LBG number density of 0.37 ± 0.21 arcmin"^ for secure LBGs and 0.56 ± 0.27 
arcmin"^ for the combined set of secure and probable LBGs. 

We modeled the photometric selection function based on the redshifted galaxy colors 
from the theoretical galaxy templates used in this survey and included photometric uncer- 
tainties modeled from the data. We combined the effects on the color selection caused by 
the characteristics of the CCDs and u'BVKl filters and assumed a uniform LBG background 
redshift distribution. The modeled data is in excellent agreement to that of the observed 
data and is well described by a Gaussian function using the parameters of the modeled data 
distribution. Moreover, a Gaussian function fits the observed redshift distributions of the 
data as well using the parameters from each of the observed distributions. We used these 
results to validate the assumption of a uniform or slowly varying LBG background number 
density from 2A < z < 3.6. The spatial information of the 211 spectroscopically identified 
LBGs from this survey will be used to measure the DLA-LBG cross-correlation with the 11 
known DLAs in Paper II. 

Following the success of the Lyman break photometric selection method demonstrated 
by Steidel et al. (1996a), several surveys for LBG candidates have extended the Lyman break 
selection method to conventional UBVRI filters (e.g.Ouchi et al. 2001; Bouche & Lowenthal 
2003; Ferreras et al. 1998; Radovich et al. 2004). These surveys rely on photometric red- 
shifts to determine galaxy spatial distributions and have a typical accuracy of Az ~ 0.2. 
The survey presented here provides deep, low resolution spectroscopy using the 10 meter 
Keck telescopes to measure the spatial distribution of LBGs at 2; ~ 3 to Az ~ 0.003. The 
339 identified spectra obtained from this survey are used to calibrate the overall efficacy of 
an LBG photometric selection method using conventional m'BVRI filters. From our spec- 
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troscopic sample and conservative identifications, we find that the criteria for the u'BVRl 
filters yield a high LBG detection efficiency (79% z > 2 identified redshifts). The inclusion of 
objects within 0.2 magnitudes of the imposed m'BVRI color criteria, made to compensate for 
photometric uncertainties, allowed for a near doubling of the number of 2; > 2 objects with 
little loss in efficiency (75% z > 2 identified redshifts). Although the use of all five filters 
with the color criteria presented in this paper may be too restrictive in detecting the full 
population of LBGs at 2; ~ 3, we find the use of a subset of these filters (m'VRI, w'BRI,m'BR) 
allow the inclusion of more objects as LBG candidates but are less efficient (69%, 63%, and 
62% z > 2 identified redshifts respectively). 

We used the spectroscopic information to offer a possible refinement to our original 
u'BVRl color criteria to increase the efficiency and provide insight into the entire population 
of LBGs at z ~ 3. Although a spectroscopic analysis of the serendipitous objects suggests 
that the original color criteria imposed by our selection method detects most of the objects 
at 2; ~ 3, they may not include all objects such as galaxies that are highly reddened. Several 
LBGs selected exclusively from a subset of the u'BVRl filters were detected outside the 
regions chosen for the original color cuts. To obtain a more complete census of objects at 
z ^ 3, modifications to the existing color cuts are suggested. 
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Observatories that made the important initial deep images in this bandpass possible and 
Kurt Adelberger and Alice Shapley for data regarding the Un filter transmission. We are 
grateful to the staffs of the Palomar and Keck Observatories for their kind assistance with 
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and Hsiao- Wen Chen for graciously allowing the use of a routine that aided in the creation of 
the stacked images. This work was partially supported by the National Science Foundation 
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Table 1. Target QSO Fields with ^ - 3 DLAs 



Field 


R.A.(J2000.0) 


Dec.(J2000.0) 


ZQSO 


'^DLA(s) 


log N(HI) 




{h m s) 


{d m s) 








LBQS0056+0125 


00 59 17.62 


+01 42 05.30 


3.16 


2.775 


21. 0^ 


PKS0336-017 


03 39 00.65 


-01 33 19.20 


3.20 


3.062 


21. 2^ 


PSS0808+5215 


08 08 49.43 


+52 15 14.90 


4.45 


2.936, 3.113 


20.9? 20.7'' 


PSS0957+3308 


09 57 44.50 


+33 08 23.00 


4.25 


3.280, 4.180 


20.5* 20.7* 


BRI1013+0035 


10 15 48.96 


+00 20 19.52 


4.41 


3.103 


21.15 


PSS1057+4555 


10 57 .56.. 39 


+45 55 51.97 


4.12 


3.050, 3.317 


20.3^''', 20.3^" 


PSS1432+3940 


14 32 24.90 


+39 40 24.00 


4.28 


3.272 


21.3* 


PC1643+4631A 


16 45 01 09 


+46 26 16.44 


3.79 


3.137 


20.78 


JVAS2344+3433 


23 44 51.25 


+34 33 48.64 


3.01 


2.908 


21.1* 



iWolfe et al. (1995) 
2prochaska et al. (2001) 

^Prom this work 
*Prochaska ct al. (2003) 
SStorrie-Lombardi & Wolfe (2000) 
SPeroux et al (2001) 
^Lu et al. (1998) 
^Schneider et al. (1991) 
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Table 2. Details of Imaging Observations 



Field 


Instrument and Date'' 


Filter 


Exp.(s) 


FWHM 


Depth AB^ 


E{B - V) 


LBQS0056+0125 


LRIS Dec 2001 




9000 


1.5 


28.9 


0.03 




LRIS/LRISb Dec 2001/Feb 2003 


B 


1530/1300 


2.0/1.2 


28.1 






LRIS Dec 2001/Dec 2002 


V 


3060/600 


1.4/1.5 


28.2 






LRIS Dec 2001/Feb 2003 


RsVR 


2800/600 


1.4/1.7 


28.3 






LRIS Dec 2001 


I 


2800 


1.5 


27.8 




PKS0336-017 


LRISb Dec 2002/Nov 2000 


u' 


7200/2400 


1.3/0.8 


29.1 


0.14 




LRIS Sep 1998 


B 


1500 


0.6 


28.9 






LRIS Dec 2002 


V 


2400 


1.7 


28.0 






LRIS Sep 1998 


R 


1200 


0.6 


28.9 






LRIS Sep 1998 


I 


1080 


0.5 


27.6 




PSS0808+5215 


LRISb Dec 2002 


u' 


4500 


1.2 


29.3 


0.04 




LRIS/LRISb Nov 2000/Dec 2002 


B 


1200/1680 


0.9/0.8 


29.2 






LRIS Dec 2002 


V 


1740 


1.2 


28.8 






LRIS Nov 2000/Dec 2002 


R 


1200/2600 


1.1/1.2 


28.8 






LRIS Dec 2002 


I 


1000 


1.1 


27.3 




PSS0957+3308 


LRISb Dec 2002 


u' 


4500 


1.0 


29.6 


0.01 




LRIS Jan 2001 


B 


3645 


1.0 


28.9 






LRIS Dec 2002 


V 


1800 


1.0 


29.1 






LRIS Jan 2001 


R 


1860 


1.0 


28.6 






Jan 2001/Dec 2002 


I 


1400/1000 


1.1/1.0 


27.8 




BRI1013+0035 


LRISb Dec 2002 


u' 


4600 


0.9 


29.7 


0.03 




LRIS Nov 2000 


B 


1200 


1.1 


28.5 






Dec 2002 


V 


1400 


0.9 


29.0 






LRIS Dec 2002 


R 


1500 


0.9 


28.9 






Lii/5 Dec 2002 


I 


1000 


0.9 


27.7 




PSS1057+4555 


LRIS Apr 2001 




9000 


0.7 


29.2 


0.01 




LRIS Jan 2001 


B 


3600 


0.8 


28.9 






LR/S Apr 2001 


V 


1800 


0.7 


28.6 






LRISJan 2001/Apr 2001 


R/Rs'= 


2320/1680 


0.8/0.7 


29.0 






LRIS Jan 2001/Apr 2001 


I 


1160/1320 


0.8/0.7 


28.1 




PSS1432+3940 


COSMIC Apr 2000 


Un 


19800 


1.0 


28.7 


0.01 




COSMIC Apr 2000 


B 


10800 


1.1 


28.8 






Lii/S Jan 2001 


V 


3200 


0.9 


28.3 






COSMIC Apr 2000 


R 


9000 


0.9 


28.5 






LR/S Jan 2001 


I 


1200 


0.8 


27.6 




PC1643+4631A 


LR/S Apr 2001 


Un 


7300 


0.9 


28.8 


0.02 




LRIS Apr 2001 


Bd 


3720 


0.9 


28.5 






LR/5 Apr 2001 


V 


1800 


0.8 


28.4 






LRIS Apr 2001 




1740 


0.7 


28.5 






LRIS Apr 2001 


I 


1200 


0.7 


28.1 




JVAS2344+3433 


COSMIC Sep 2000 




32400 


1.4 


29.2 


0.08 




COSMIC Sep 2000 


B 


8400 


1.2 


28.8 






Li?/5 Jan 2001 


V 


2400 


1.2 


28.9 






COSMIC Sep 2000 


R 


7980 


1.0 


28.7 






Lii/5 Jan 2001 


I 


1200 


1.3 


26.7 





''LRIS: The Low Resolution Imaging Spectrometer; LRISb: LRIS witii the new blue sensitive CCD chips; COSMIC: 
The Carnegie Observatories Spectrograph and Multi-Object Imaging Camera. 

^la sky fluctuations per seeing element. 

''The Rs filter was used instead of, or in addition to, the Rjc in this field. 

''Johnson B filter convolved with the transmission of the 460nm dichroic installed on LRIS 
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Table 3. LBG Color Candidates 



Field 


Area'' 


«'BR(I) 


u'VRI 


u'BVRI 


u'YRI LBG 


u'BVRI LBG 




(circmin^) 


Candidates 


Candidates 


Candidates 


Density'' 


Density"^ 


LBQS0056+0125 


41.96 




77 




1.84±0.2ld 


d 


PKS0336-017 


40.95 


87 








e 


PSS0808+5215 


45.01 


61 


54 


42 


1.20±0.16 


0.93±0.14 


PSS0957+3308 


39.75 


47 


34 


29 


0.86±0.15 


0.73±0.14 


BRI1013+0035 


37.21 


45 


19 


16 


0.51±0.12 


0.43±0.11 


PSS1057+4555 


41.49 


68 


81 


46 


1.95±0.22 


1.11±0.16 


PSS1432+3940 


87.16/57.40^ 


165 


79 


42 


1.38±0.15 


0.73±0.11 


PC 1643+4631 A 


38.08 




91 




2.39±0.25'' 


d 


JVAS2344+3433 


93.14/59.62* 


155 


51 


40 


0.86±0.12S 


0.67±0.118 



^'Common area to all filters. 

''Number of LBG candidates per arcmin^ iising u'VRI filters for compajrison to previous surveys. 
"^Number of LBG candidates per arcmin^ using u'BVRI filters. 

"^LBG candidates were chosen using u'VRI filters. 
''LBG candidates were chosen using u'BRI filters. 

*Area common to the 3 COSMIC filters (u',B,R) and area common to all 5 filters [COSMIC (u',B,R) and 
LRIS (V,I)], respectively 

SLBG candidates were chosen using ti'BR filters. 
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Table 4. Details of Spectroscopic Observations 



Field 


Date 


Slitmask" 


Area** 


Exp.(s) 




FWHMd 


LBQS0056+0125 


Nov 2003 


q0056I15 


41.3 


4500 


37 


0.8 




Nov 2003 


q0056rl5 


41.3 


2400 


37 


0.9 


PKS0336-017 


Feb 2003 


295n0336 


41.3 


4500 


41 


1.1 




Feb 2004 


405q0336 


41.3 


3600 


28 


0.8 




Feb 2004 


420q0336 


41.3 


3600 


27 


0.9 


PSS0808+5215 


Feb 2003 


q0808nl5 


41.3 


5234 


27 


0.9 




Nov 2003 


122q0808 


41.3 


4200 


30 


0.9 




Feb 2004 


q0808cl5 


41.3 


4500 


25 


0.8 


PSS0957+3308 


Feb 2003 


180p0957 


41.3 


4500 


21 


0.9 




Nov 2003 


100q0957 


41.3 


3600 


25 


1.0 




Feb 2004 


100r0957 


41.3 


4500 


19 


0.9 


BRI1013+0035 


Feb 2003 


ql013nl3 


41.3 


4500 


26 


0.9 




Nov 2003 


ql013rl5 


41.3 


3600 


22 


1.0 




Feb 2004 


ql013rl5 


41.3 


4500 


22 


1.0 


PSS1057+4555 


Dec 2002 


rl057nl0 


41.3 


3600 


36 


1.0 




Dec 2002 


11057tl0 


41.3 


6000 


33 


1.1 


PSS1432+3940 


Feb 2003 


ql432pl5 


41.3 


4500 


31 


1.2 


PC1643+4631A 


May 2003 


ql643rl5 


41.3 


4800 


30 


0.7 




Feb 2004 


245ql643 


41.3 


3400 


40 


0.8 


JVAS2344+3433 


Nov 2000 


q2342110 


31.5 


7200 


18 


0.8 




Nov 2000 


q2342cl0 


31.5 


7200 


23 


0.9 




Oct 2001 


q2342rl5'= 


31.5 


7200 


21 


1.2 




Dec 2002 


2:-!12i)2S0 


11.:-! 


7200 


:-!0 


1.5 



'^Name of slitmask for reference. 

''Approximate area covered by the slitmasks in arcminute^. Slitmasks overlapped in 
most fields. 

"^Number of slitlets milled to obtain object spectra. Some objects were observed on 
more than one mask to obtain a higher signal-to-noise ratio. 

Seeing FWHM taken from alignment stars in arcseconds. 

''These images were degraded due to an optical oil leak onto the collimator of LRIS. 
Objects from this slitmask were not included in the sample. 
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Table 5. Spectroscopic Identifications 



Field 


No6." 


Starts 




Low-Zprob'^ 


Kigh-Zprob" 


LBQf 


QSQs 


Unidentified'' 


(Category) 




(0) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


LBQS0056+0125 


73 


3 


12 


26 


9 


12 


1 


10 


PKS0336-017 


73 








5 


8 


14 


1 


45 


PSS0808+5215 


59 





4 


11 


13 


14 


1 


16 


PSS0957+3308 


48 


4 


6 


6 


5 


16 


1 


10 


BRI1013+0035 


39 


3 





5 


7 


16 





8 


PSS1057+4555 


63 





2 


2 


5 


20 


1 


33 


PSS1432+3940 


31 





5 


5 


5 


2 





14 


PC 1643+4631 A 


73 





3 


6 


13 


29 





22 


JVAS2344+3433 


70 


1 


15 


4 


9 


9 





32 


Totals 


529 


11 


47 


70 


74 


132 


5 


190 



''Total number of LEG candidates observed spectroscopically 

''Number of identified Galactic stars 

"^Number of spectroscopically secure z <2 galaxies 

''Number of objects with probable z < 2 identifications 

''Number of objects with probable z > 2 identifications 
^Number of spectroscopically secure 2 > 2 galaxies 
^Number of identified QSOs with z>2 
''Number of unidentified objects 
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Table 6. Color Selection Criteria Spectroscopic Efficiency 



Filters'- 


Star^ 


LoW-Zga('= 


Low-Zprob'^ 


High-Zprob" 


LBG* 


QSOS 


Unidentified'' 


z > 2' 


u'BVRI 





3 


10 


12 


35 


1 


15 


0.79 


u'BVRI ±0.2j 


5 


9 


16 


24 


61 


3 


41 


0.75 


u'VRI 


1 


3 


5 


8 


11 


1 


15 


0.69 


w'VRI ±0.2j 


3 


18 


27 


18 


37 


1 


47 


0.54 


w'BRI 





4 


12 


12 


14 


1 


45 


0.63 


u'BRI ±0.2j 


2 


5 


21 


19 


24 


1 


69 


0.61 


u'BR 





6 


4 


8 


8 





14 


0.62 


u'BR ±0.2j 


1 


15 


4 


9 


9 





28 


0.47 



'^Objects meet the criteria of the listed filters only 

''Number of identified Galactic stars 

"^Number of spectroscopically secure 2 < 2 galaxies 

'^Number of objects with probable z < 2 identifications 

''Number of objects with probable z > 2 identifications 

'^Number of spectroscopically secure z > 2 galaxies 

SNumbcr of identified QSOs with z>2 

''Number of unidentified objects 

'Fraction of secure and probable objects with z > 2 

JColor criteria extended 0.2 magnitudes to compensate for photometric uncertainties 
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10000 



Fig. 1. — Transmissions of the filters used in this survey. The dotted curves indicate the 
effective transmissions corrected for atmosphere attenuation and the appropriate CCD quan- 
tum efficiencies. The [/„ and filters were used in certain fields and are denoted by the 
dashed curves. 
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Fig. 2. — Number counts versus magnitude in five fields with u'BVRI imaging. This method is 
used as a check for photometric completeness. The upper histogram shows the number of source 
detections versus R magnitude (hatched region) binned in AR= 0.2. Other filters are plotted versus 
their respective magnitude with line styles indicated in the legend. The spectroscopic magnitude 
limit of LRIS is denoted by a vertical dotted line at R= 25.5. The central plot displays a logarithmic 
fit (solid line) to the number counts versus R magnitude. The lower plot displays the logarithmic 
fit to all five filters for direct comparison. In this plot, the error bars have been removed for clarity. 
The slope of the logarithmic fit flattens for consecutively redder filters. 
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Wavelength A 

Fig. 3. — Graphical depiction of the Lyman break photometric selection. The filter transmissions 
of the three filters (u'BV) are shown with their transmissions before (solid curves) and after the 
correction for CCD quantum efficiency and atmospheric attenuation (dotted curves). Vertical 
dashed lines indicate the effective wavelengths of each filter. Overlayed is a spectrum of a bright 
LBG from the data with its wavelength array shifted to represent an arbitrary LBG at z = 3.0. 
Short-ward of 912A in this spectrum is shown as zero flux. The placement of the break in the 
continuum of the galaxy at the rest-frame Lyman limit (9 12 A) and at rest-frame Lyman a can 
be seen with respect to the filter responses. This particular LBG displays a Lyman a emission 
spike (near 4900 A). Our color criteria was chosen to detect a distribution of objects that have the 
Lyman a break centered at z = 3.0 causing their Lyman limit break to be redshifted to within the 
passband of the u' filter. 



- 36 - 



I 




> 

I 








B-R 



V-R 



> 

I 




2 3 



4 



V- 



Fig. 4. — Color-color plots showing the expected colors of a starbursting galaxy (dashed green 
curve) and an elliptical galaxy (dotted blue curve) as they are evolved from z = 0. The star- 
bursting galaxy is indicated by two solid red curves from 2.5 < z < 3.5. The separate red curves 
indicate decrements in the continuum of the galaxy caused by a reasonable range of Lyman a forest 
absorption. Star symbols illustrate the colors of Galactic stars from the Gunn and Stryker (1983) 
catalog. The LBG candidate selection regions are the boxed areas in the upper-central or upper-left 
portions of the plots. Low redshift reddened elliptical galaxies and low-mass Galactic stars are the 
main source of interlopers. 
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Fig. 5. — Color-color plots of the data from the fields with u'BVRI color selection. All 20.0 < 
R < 25.5 objects are denoted by red points. Those objects meeting all n'BVRI color criteria are 
indicated by green diamonds and reside in the boxed selection regions in the upper central, or upper 
left portions of the plots. Black asterisks denote objects meeting all u'BYHI color criteria and have 
secure z > 2 redshifts. The boundary regions surrounding the selection regions reflect a relaxation 
of the color selection criteria by 0.2 magnitudes to compensate for photometric errors. Most of the 
objects selected are not detected in the n'-band, therefore the values shown are the 2a lower limit 
to their (u' — B) and (n' — V) colors. 
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Fig. 6. — Composite LBG spectra shifted to the rest-frame. Each spectrum consists of multiple 
z > 2 galaxy spectra that exhibit similar properties. Expected interstellar absorption and emission 
lines are labeled and indicated by the dotted lines. Top: Composite spectrum of 25 LBG galaxies 
that display prominent Lyman a in emission. Interstellar absorption lines are typically weaker 
in this population. A detectable break in their continua exists short-ward of 1215A. Bottom: 
Composite of 15 LBG spectra displaying Lyman a in absorption. These objects display stronger 
interstellar absorption lines and a slightly reddened continuum. Fluxing beyond rest-frame 1500A in 
each composite spectrum is less reliable for reasons discussed in the text. 
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Fig. 7. — Histogram of the R magnitude distribution of the spectroscopically identified objects 
with z > 1. The secure LBGs (hatched region) are shown with probable LBGs (sohd region) 
overlayed. 



-40- 



CO 

O 

m 



03 

E 



m 
O 

m 

I 

o 
O 



0) 



o 
m 



0) 

E 



0) 

c 

JD 
E 

o 
o 



14 

1 2 
10 
8 
6 
4 
2 



15 



10 



Lya absorption — hatched 
Lya emission - solid 



22 




23 24 
R magnitude 



Lya absorption — hatched 
Lya emission — solid 



22 



23 24 
R magnitude 



Lya absorption - hatched 
Lya emission - solid 





23 24 
R magnitude 

Fig. 8. — Histograms of the R magnitude distributions of the emission-identified (solid region) and 
absorption-identified (hatched region) LBGs. The upper and central histograms plot the secure 
and probable LBG distributions respectively. The lower histogram plots the distributions of the 
combined set of LBGs. 
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Fig. 9. — Redshift distributions of LBGs displaying Lyman a in absorption (hatched region) and 
emission (sohd region). The secure LBG identifications are presented in the upper histogram and 
probable LBG identifications are presented in the center histogram. The lower histogram displays 
both secure and probable LBGs combined. The broad similarities in the absorption-identified 
redshift distributions reinforces the overall probable redshift assignments. 
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Fig. 10. — Real and simulated redshift histograms. The top histogram represents the redshift 
distribution of 137 LBGs from the data with secure spectroscopic redshift identifications. The 
central histogram includes an additional 74 objects with probable redshift identifications. The 
LBG selection function is depicted by the solid curve in both cases. The lower histogram displays 
the results of a simulation of the observed redshift distribution accounting for the characteristics 
of the instrument, filters, and observed photometric uncertainties. The simulation assumes a flat 
background distribution of LBGs over the redshift path 2.4 < z < 3.6. The normalized residuals to 
the Gaussian function fits with respect to the expected Poisson fluctuations per bin are indicated 
by dashed curves. 
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Fig. 11. — Data from five fields with ii'BVRI color selection are plotted to re-assess the original 
photometric selection criteria (dashed lines). The z > 2 LBGs meeting u'BVRI color criteria are 
marked with black diamonds. Those meeting only u'BRI or u'VRI color criteria are marked with 
blue squares or green triangles respectively. Interloping objects with z < 2 redshifts from the data 
meeting m'BVRI color criteria are marked with red asterisks whereas those meeting only u'BRI or 
u'VRI color criteria are marked with red plus '+' signs or 'X's respectively. For the objects that are 
non-detections in the u'-band, the points represent the 2a lower limit to their (u' — B) and {u' — V) 
colors. Spectroscopic analysis of confirmed LBGs shows that a fraction of the galaxy population 
at z ~ 3 is missed using the full u'BVRI criteria. To probe a more complete sample of galaxies at 
z ~ 3 at the risk of lower efficiency, a possible refinement to the u'BVRI photometric selection is 
shown (solid lines) in the upper-left portions of the plots. 



